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The existence of superfluidity in a 3D Bose-gas can depend on boundary interactions with channel 
walls. We study a simple model where the dilute moving Bose-gas interacts with the walls via 
hard-core repulsion. Special boundary excitations are introduced, and their excitation spectrum is 
calculated within a semiclassical approximation. It turns out that the state of the moving Bose-gas 
is unstable with respect to the creation of these boundary excitations in the system gas + walls, i.e. 
the critical velocity vanishes in the semiclassical (Bogoliubov) approximation. We discuss how a 
condensate wave function, the boundary excitation spectrum and, hence, the value of the critical 
velocity can change in more realistic models, in which "smooth" attractive interaction between the 
gas and walls is taken into account. Such a surface mode could exist in "soft matter" containers 
with flexible walls. 
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I. INTRODUCTION 

Recent experiments on Bose-Einstein Condensation (BEC) in magnetically trapped gases 0, and excitons in 
semiconductor crystals and nanostructures ||, || have made the subject of BEC more vital, more interdisciplinary. 
Many new questions appeared naturally as the understanding of the process of BEC progressed ||. However, some 
"old" problems - such as the kinetics of BEC, the nature of superfluidity, the critical velocity problem, etc. - still 
remain the subject under consideration R], (especially for new physical objects where the state of BEC was recently 
demonstrated Q). Although the critical velocities are one of the first difficulties to be encountered in the study of 
superfluidity, they are still the least understood aspect in the theory [||, ||, poj . 

This article is motivated mainly by the problem of critical velocity(-ies) in exciton superfluidity. It has been found 



experimentally ||], Jll|, |12| that a cloud of condensed excitons moves through a crystal with some constant velocity 
and some characteristic shape of the density profile. Several theoretical explanations of this anomalous transport have 
been put forward |13| ], p[ , J15[ ] . In spite of the fact that these explanations are based on different assumptions, there 
are several common ideas in the background of all these theories. For instance, it is the notion that interaction with 
a lattice is very important (if not to say crucial) in the BEC of excitons [^5) , Q , jlTj . 

However, there are many outstanding questions that remain the subject of discussion [ fLs] , One of such questions, 
for example, is the superfluid nature of exciton anomalous transport. In fact, it is not clear how the exciton condensate 
"feels" the boundary of the crystal (via interaction with surface phonons, e.g.,) or the impurities and other lattice 
imperfections that can bound an exciton. Generally speaking, clarification of the role of these friction sources may be 
essential for the understanding of the exciton superfluidity and the nature of critical velocities. 

We approach the critical velocity problem by working out a simple model in which dilute 3D Bose-gas moves in a 
channel and interacts with the walls of this channel. The walls are modeled as two 3D solid bodies with well-defined 
boundaries. Although we take into account repulsive interaction between the particles of the gas, the proposed model 
cannot describe, for example, the superfluid He, which is a Bose liquid with strong interparticle interaction. Yet 
this is not the aim of this article. The main goal of this study is to explore the space of manaeuvre appearing in 
the framework of the well known simple models, such as the weakly nonideal Bose-gas, if we switch on the gas-wall 
boundary interaction. 

We show that the existence of the repulsive interactions between the Bose-gas and the channel walls leads to the 
essential reduction of the critical velocity of the superflow. The finiteness of the Landau critical velocity for the bulk 
(i.e. Bose-gas) excitations turns out not to be a sufficient condition of superfluidity. Note that in the present work 
we investigate the boundary excitations. Although the breakdown of the superfluidity is assumed to be accompanied 
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by vortex emission, we leave for future studies the questions of the vortex formation and their dynamics in the case 
when interaction with walls is taken into account. 



II. CRITICAL VELOCITY PROBLEM 



A closed system cannot undergo an inner macroscopic motion in thermodynamic equilibrium. Once such a motion 
is present, the system must evolve toward an equilibrium state. However, unless this transition is kinematically 
prohibited, (i.e. incompatible with conservation laws), the macroscopic motion is sustained. Such is the case with a 
small object moving without any viscous drag in stationary superfluid fl9|| , f20f . This object is assumed to have no 
inner degrees of freedom, so that its momentum and energy depend only on the velocity v of the center of mass. 

If the conservation laws for the creation of an excitation with the energy e gas (fc) and the momentum k in a fluid 
(gas) lead to |§- || 

Cgas(fc) > 0, in the object reference frame, (1) 

the particle will continue to move without any experience of drag forces. Condition ([!]) is known as the Landau 
criterion of superfluidity. In fact, Eq. ([!]) holds that 

v < = (e gas (fc)/fc) min . (2) 

Here Vt, is the Landau critical velocity. 

Formula (||) is in agreement with experiments performed with the semimicroscopic objects moving in the liquid 
helium 

Formula ([!]), (taken in the channel reference frame), is employed as a criterion of Bose-gas superfluid flow in 
channels. In that case, the channel walls are regarded as a massive macroscopic body in the above consideration (i.e. 
the walls act as some source of perturbations on the gas flow), and the final result is formulated in the form (||). 

On experiments with liquid helium flow, however, the registered values of the critical velocities turn out to be much 



smaller than «l- Moreover, the critical velocity depends on the channel dimensions 22 . The fact that the liquid 
superfluid helium could not be treated as a dilute Bose-gas is believed to be the main reason for this discrepancy. It is 
generally assumed that the critical velocities are related to the appearance of quantized vortex lines in the superfluid. 
The Landau criterion (Q) applied to the vortex excitations |23| l can explain the critical effects in circular geometries. 
However, it cannot account for the drastically different critical velocities for rotation and linear flows p4j . 

A superflow of a dilute Bose gas, described by nonlinear Schrodinger equation [p5| , has been studied recently in 
different geometries with the use of direct numerical methods ^6|- [ p7| . It has been observed that the distinct critical 
velocity is linked to the emission of vortices. This velocity turns out to be equal to the Bogoliubov |2£| velocity of 
sound propagation in the gas. It is in good agreement with the criterion (^) since the Dirichlet boundary conditions 
(for Bose-gas wave function) are imposed on the channel walls ^7j. This means that the walls, being considered as 
rigid immovable bodies, have no degrees of freedom. As a consequence, the arguments leading to formulas (0),(§) can 
be used. 

In reality, the channel walls have a large number of degrees of freedom. Indeed, short- and long-range forces between 
a particle and a surface, boundary and interface phonons are well known subjects in the surface physics p9| , | B(j |. This 
means that the (superfluid) Bose-gas is coupled with the channel, in which the gas moves. Then special boundary 
excitations can exist in the system of Bose-gas + channel walls because of the coupling between, say, the surface 
phonons of the walls and the Bogoliubov phonons of the Bose-gas. Therefore, the Landau criterion in the form (|l|) 
cannot be applied; it has to be modified. To get an analog of it we use the laws of conservation, taking the (inner) 
walls' degrees of freedom into consideration. The superflow can exist, provided the following condition holds: 

e(fc) > 0, in the channel (laboratory) reference frame, (3) 

where e(k) is the energy of any elementary excitation of the whole (gas + walls) system. 

Condition (§ means that the state of the system can not be changed, since the occurrence of any number of 
elementary excitations leads to the increase of a total energy, but the latter is prohibited by the low of conservation. 
The excess of the momentum is "absorbed" by the motion of the center of mass of the walls. The energy is not 
actually changed by this motion (in the channel reference frame) because of a large mass of the walls and their zero 
initial velocity. 
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III. BOSE-GAS WITH BOUNDARY INTERACTIONS 



We study the model of the dilute Bose-gas in the channel of the width 21 (in y-direction) and of infinite length (in 
x- and z-directions). The channel walls occupy \y\ > I part of the space (see Fig. 1). The general structure of the 
Hamiltonian is the following: 



h = ff gas (V>, ^) + ffph^gi, 7Ti) + fli n ti(gi, V> f V>)- 

+ -ffph 2 (<?2, K 2 ) + H iat2 (q2, fti>), 



(4) 



where ip is the Bose-gas field operator, q is the displacement field operator of a wall, tt is the momentum density 
operator conjugate to q, and the indexes 1,2 correspond to the upper and lower part of the channel respectively. In 
the model being considered, the Bose-gas Hamiltonian has the following form: 

-ffgas = Jft(r) (-5^) $(r) dr + J ^S(r - r')ft (r)ft (r')^(r')^(r) drdr 1 , 
where v > is the interparticle interaction constant |20| , and the wall Hamiltonian can be written as follows 



H 



P h 



2p 



+ djq k (r)\ jkln diq n (r) dr, 



where the tensor Xjkin describes the elastic properties of the channel walls. 

We derive the excitation spectrum using the techniques of semiclassical approximation (cf. 28 1, |||], p3|). 
Expanding the field operators near certain classical solutions, i.e. ip = ipo + Sip and q — qo + Sq, we represent the 
Hamiltonian in the form 



H = H + UH 2 



(5) 



where Ho stands for the classical part of H. Note that ipo ^ indicates existence of a condensate in the moving 
Bose-gas, whereas go ^ appears in this model mainly to satisfy boundary conditions (see below). The Hamiltonian 
H 2 in (||) is bilinear with respect to the field operators. As a consequence this (semiclassical) Hamiltonian can be 
reduced to the normal form 



Hi 



Ml* 



const, 



h,,b] 



= 5,. 



IJ1 



b i} bj 



= 0. 



(6) 



The quantum Hcisenberg and the classical Poisson-Hamilton equations of motion for the field operators (functions) 



i(dh/dt) = 



HA 



(dh/dt) = {H, h} , 



(7) 



have the same form, if in (j^) we neglect the terms of power in Ti greater than one. These equations are linear with 
respect to the field variables. 

It follows from (^), (0) that the excitation energies u>i are equal to the characteristic frequencies of these equations. 
Thus, to determine the semiclassical energy spectrum of the system we need to find the characteristic frequencies of 
the classical field equations for Sip and Sq linearized around a proper stationary solution tpo, qo. (This means that 
Sip and Sq, originally of the operator nature, can be treated as a c- number and a real number function respectively). 
Then the gas state can be characterized by the Bose wave function ip(r, t) — tp(x, y, t), while the state of the walls is 
determined by the displacement field q(r,<) = (q x (x,y,t), q y (x,y,t), 0). For simplicity, the system is asumed to be 
homogeneous in the z-direction. 
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FIG. 1. The Bose gas moves with the velocity v in a channel of the width 21. The profile of the stationary wave function 
4>{y) = IV'oC*, V, t)\ is shown in the center, where A is the coherence length. The top and bottom curves depict the boundaries 
of the walls interacting with the gas, while the horizontal bold dashed lines correspond to the unperturbed boundaries, and 
q y (x,±l,t) are deviations from this equilibrium. 



The interaction between the wall and gas atoms is given in our model by the sum of hard-core repulsion and some 
"smooth" potential U = W(r, q). The hard-core repulsion makes the walls impenetrable for the gas. It follows 

that the wave function of the Bose-gas vanishes along the actual boundaries y± = ±1 + q y (x,±l,t) (see Fig. 1). 

ip(x,±l + q y (x,±l,t),t) = 0. (8) 

The other boundary conditions, namely 

h 2 

(Jyy(x, ±1, t) = —9yljj* dy4>(x, ?/ = ±1 , t) , (9) 

a xy (x,±l,t) = 0, (10) 

correspond to the equality of the forces between the solid and gas on the boundary. In (||) , (|l(]) 0y denotes the wall 
stress tensor and to is a mass of the gas atom. The pressure of the gas (r.h.s. of @) is equal to the normal component 
of the stress tensor (l.h.s. of (g)). Eq. (|l^) follows from the fact that the tangent stress vanishes on the boundary. 
The wave function ip of the repulsive Bose-gas satisfies the nonlinear Schrodinger equation p5|] 

ihd t + —A-vi)*ip\i> = V(r,d j q k )i>, \y\ < I, (11) 

while the wall dynamics obeys the hyperbolic equation for the displacement field q J34j, 

pd 2 qi = djCTij - diW(r, ^*V0, |vl > h (12) 

(p denotes the wall mass density). The potentials V and W depend on the "smooth" part of the gas- wall interaction 
U; they must vanish if U = 0. To make our model as simple as possible, we set hi — a priori. Then, dynamics 
is described by eq uations (pd|), with constant coefficients, whereas the hard-core interactions fix the boundary 
conditions (||)-(10). Note that these conditions imply that the walls have a finite compressibility, K = —V(dvp)s < °°- 
Therefore, even though U — 0, the wall and Bose-gas excitations can be coupled because (d y ip) q y (x 1 ±l) in (||) and 
Oyy ~ Kd y q y (x, ±Z) in (|^) are finite and time dependent. 

We suppose that the walls are isotropic (this assumption does not affect results qualitatively, simplifying our 
calculations), so that @ 

O-y = P c 2 (d jqi + d iQj + {fi- l)%(Vq)) , (5 = (c 2 /c 2 ) - 1, (13) 

where q , Ct are longitudinal and transversal sound velocities respectively. It is convenient to rescale variables in such 
a way that the spatial coordinates are expressed via the coherence length A units while the flow velocity is measured 
in terms of the Bogoliubov sound velocity cb J^g] : 

X = H/ ytlSPgas, C B = ytvpgas/m. 
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Here, p ga , s stands for the bulk gas density. Then Eqs. (|1 l|),(p_2[) become 

(id t + (1/2) A - r V) fl>(x, y, t) = 0, (14) 

{c t ~ 2 9 t 2 - A - /3V(V-)} q(s, y, t) = 0, (15) 

where the sound velocity, time and q are measured in the units of cb, A/cb and A respectively. 

In the stationary regime, the Bose gas moves uniformly in the rr-direction with the velocity v. As the system is 
homogeneous in the x and z directions, |i/>o( r )i)| = (j>{y) and qo(r,i) = qo(y)- The walls are deformed only in the 
y-direction, since the tangent stress vanishes in the stationary regime qo(j/) = (0, qo y {y), 0). The corresponding 
solution of (O), ( p"5| ) is given by 

M*, y, t) = <P(y)e ivx e~ int , 0(±O = o, 0'(±z) ^ o, (16) 

q 0x = 0, q 0y = Q(y) = ±const(y ± I), (y -» ±Z), Q(±Z) = 0, Q'(±l) £ (17) 

It follows from (fl6|), ( |l4| ) that f2 = /2 + tj-, ft = 1 (this corresponds to the value of a chemical potential /i = vp sas /m 
at T = 0), and 0(z/) satisfies the following equation fl35| : 

-^" + 3 = </>, <Kv) = <K~v). (18) 

The parity of in ( fl8| ) and the boundary conditions in ( p"6| ) are obtained from (p|)-(^0|) (see Fig. 1). 
We follow the procedure (H)-(Q) expanding the field variables around the stationary solution ([if]), ( p7| ) 

ip = (<Kv) + £(* - «t, y, t) )e ira e"*^/ 2 )*, 
= (x(x - vt,y,t), q y = Q(y) + ( y (x - vt,y,t). 

Substituting these expansions into (|I^,(|l5|) we get the following linear differential equations for fluctuations 

^ • — ~ " s 2 1 t Al„.\1t* 



id t + 2 A + 1 - 24>{yYj £ - <t>{yYe = 0, \y\ < I, (19) 
{c t - 2 (d t -^) 2 -A-/3V(V-)}C = 0, \ y \>l, (20) 

€ = Z(x,v,t), C = C(x,y,t), 

written in the reference frame moving with the Bose-gas, x — > x' — x — vt. 

One of the advantages of setting U — is the possibility of using the exact solution of Eq. (Il8|) with the boundary 



conditions (16) 



<f>(y) = sn(y + I, q), 

where sn(y, q) is the elliptic sine f36[| , the parameter q is chosen to fit the boundary conditions and the following 
condition holds sn(y, q) — > tanh(y) if I — > oo. (Notice that the dimensional condensate wave function can be written 
in the form <f>d(y) = const ^p gas /msn( {y + l)/X, q) ). 

Nontrivial excitations can not propagate over the wall region far from the boundary. Indeed, the equation ( ]20| ) 
has the constant coefficients and, hence, the dispersion law for such excitations coincides with the phonon one (i.e. 
corresponding asymptotical solutions describe propagation of the ordinary sound waves far from the boundaries). 
Therefore, we have to look for the solution of ( po|) decreasing in y — > ±oo directions. We assume that the experimen- 
tally discovered dependence of the critical velocity on the canal width | p2[ , v c ~ l~ n , n ~ 2, is a hint to search for the 
special type of excitations, in which two boundaries of the canal can contribute coherently. 

Such a solution of (|l9|), ( pp| ) can be written in the form 

i = Xi(y)sin(kx-u)t)+ix2(y)cos(kx-uit), \y\ < I, (21) 
(x = n(y) cos(kx - cut), Q y ^ r 2 (y) sm(kx - ut), \y\ > I (22) 



with 



i(y) = A, exp(— #s|y|) + B t exp(-r)\y\). (23) 
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Two exponential terms in (|23|) correspond to the different polarizations of the boundary excitations. The char- 
acteristic values k = \/k 2 — oj 2 /c 2 , r/ = ^Jk 2 — uj 2 /cf are eigenvalues of the ordinary linear equations obtained by 
substitution of ( p2[ ) into (ppf). Note that the ansatz (|2l|), p^) is equivalent to the Bogoliubov u-v transformation 
generalized to a nonuniform case pH , | 02| , |33| and coupling with the surface phonons: 

e iflt 6^{x, y, t) = u k {y) e^-'-W*) + e i(-fcx+u.(*)t) 
5q(a ; ,y I t) = C fe (y) e^""^ + c.c. 



Then the operators &t, (6^) that create (annihilate) the boundary excitations (21-23|) in the diagonalized Hamiltonian 
(^) can be represented by the linear combinations of the Bose-gas field operators, 5ifi and Sffi, and the displacement 
field operators, Sq and dir. 

We linearize the boundary conditions (|s|)-(|lO|) according to the method of the semiclassical approximation (| ) , 
neglecting the terms of power greater than one in £, Together with the proper solution (|2l|),(^2|) of Eqs. (|l9|), (EC), 
conditions (p|)-(^0[) determine values A = A(k,u>), B = B(k,uj) in ( p3| ) and the boundary conditions for Xx,2'- 

—) =Tk-y(z), X2(±0 = 0, z = r%i, (24) 



pc 2 /p gas c 2 B 1_ ( 4VT — z _ (2-zf 



> i. 



The value e(k) in (M) 



e(k) = uj(k) — kv 



equals the boundary excitation energy in the channel reference system (g). 

It follows from ( pl| ) and (|l^) that the variables Xi,2(y) an d u> satisfy the following system: 

L lX x{v) + 2io X 2(y) = 0, Lx = d 2 - k 2 + 2 - 6<P(y) 2 , 

L 2 X 2 (y) + 2 WXl (y)=0, L 2 = d 2 -k 2 + 2~20(y) 2 . (25) 

Note that in view of the parity of tfi(y) (see (|l8|)) and the boundary conditions (|2~i|), solutions of ( p5| ) can be either 
symmetric or antisymmetric with respect to y. 

According to the criterion (^), a breakdown of superfluidity occurs at such a value v if there exists such a k ^= 
that e(fc) = 0. Then the argument z of 7 in (|24|) vanishes and the boundary conditions do not depend explicitly on 
u>, i.e. {x'x I Xx)\y=±i = ^^7(0) = Tconst. In principle, this fact makes it possible to calculate v c without finding any 
final expression of uo(k). Indeed, one has to solve Eqs. ftlEj) with uj = kv c and fixed boundary conditions. 

In this simple model, however, it is possible to calculate the dispersion relation uj(k), at least in the k — * limit 
|57j. We look for the symmetric solution of (|25|), expanding Xx,2(y) in powers of uj 2 : 

(0) 2 (1) , 1 (o) . 2 (1) , \ 

Xx = Xx + & Xx + • • • , X2= ^(x 2 + w x 2 + ■ ■ •)■ 

The functions X ^\ symmetric in y satisfy the following recurrence relations 

L^ 0) =0, x ( x°H±l) = l, 

L^+xt 1 '* =0, Xi i) (±0=0,i = l,2,... 

L2xf + xf = 0, X ( 2 } (±1) = 0, » = 0, 1, 2, . . . . (26) 

The analytic study of (^6|) seems to be difficult. Instead we proceed numerically psj ; we obtain the eigenvalues 
w = oj(fc) by solving (|2^) recursively and imposing the boundary conditions (|24|) on x- 
The result reads 



= oV(7 + <5)fc 3 , as fc -> 0, (27) 

where a > and 5 > are some bounded functions of Z, 7 = 7(0) ~ pc 2 / p gas c B . The dependence of log(w(fc)) on log(/c) 
is shown on Fig. 2. Note that only the solutions with finite 7 > 1 have the physical meaning. For "conventional" 
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systems, such as a dilute Bose gas inside a solid container, the value of 7(0) is very large, 7(0) ~ 10 7 ~ 8 or even more. 
Moreover, the validity of the long- wavelength / low-energy approximation implies kj < 1 and and the relevant 

wavelengths are unphysically huge. However, for the "soft matter" substances with flexible walls, 7 can be of the 
order of 10 2 ~ 3 and the distance between the walls can be 21 > A. Then, beginning from the wavelengths of the order 
of (10 25 ~ 3 )A, we are within the 'fc — + 0' limit and u>(k) < 10~ 3 /2. (On Fig. 2, we present also the curves with 7 < 1 
because of similarity between our result and the dispersion relation of capillary waves on the interface between liquid 
and gaseous He, the so-called "ripplons" If the channel walls were rigid and incompressible, that corresponds 

to 7 = 00 in (|27|), the inhomogeneous surface excitations introduced in this study just do not exist as a well-defined 
object. 

It is easy to conclude that the semiclassical critical velocity is zero in this model, since for any v there exist k 7^ 
such that e(k) = u)(k) — kv < 0. This means that the model, in which all the gas -boundary interactions are reduced 
to the hard-core repulsion, predicts (in the semiclassical (Bogoliubov) approximation) an instability of the Bose-gas 
current state in relation to occurrence of boundary excitations. However, whether the damping of superflow can 
happen via the energy transfer from the 3D condensate to the boundary localized modes is an open question, which 
cannot be answered in the framework of models with the supcrfluid density p s — p gas , T = 0. More sophisticated 
models of superfluidity, such as the two-fluid hydrodynamics p2j , should be used to describe the kinetics of damping. 
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FIG. 2. The low-momentum spectrum u)(k) of boundary excitations. The canal width I is measured in the coherence length 
units; 7 ~ pc 2 t / p sa _ s c\. 

IV. DISCUSSION 

In our opinion, the benefit from finding such a dispersion relation (see Fig. 2) in the framework of the proposed 
oversimplified model is the possibility of concluding that the boundary interactions can play a key role in the super- 
fluidity, namely, by reducing substantially the critical velocity of the superflow. If that is the case, it is reasonable to 
generalize our model to make it more realistic. 

Recall that we have neglected the "smooth" part U of the gas- wall interactions. The models with repulsive U seems 
to be qualitatively similar to the one under consideration. We believe that they also lead to the zero semiclassical 
critical velocity, since nothing can apparently prevent the energy transfer from the gas flow to the walls (in this 
approximation). It is important to note that even if the higher quantum corrections yield nonzero critical velocity, 
the latter will be of much smaller value than the (semiclassical) Bogoliubov-Landau velocity found by the numerical 
simulations based on the nonlinear Schrodinger equation p7| . 

The situation might be different if the "smooth" part U of the interaction between the gas and wall atoms were 
attractive. For example, we can take the Hamiltonian of gas-wall interaction in the Deformation Potential form 

H int - / a(r, r') Vq(r) ^(r') rfrrfr', (28) 
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where r and r' change in the wall area (\y\ > I) and in the channel area (\y'\ < I) respectively, and the function 
er(r, r') describes atom-lattice interaction. (At least two new parameters, which control the "smooth" part of gas-wall 
interaction, have to appear in the model with W ^ 0: one, some characteristic value of energy, and, two, a length 
scale). Then the equations for the classical parts of ip and q, namely (p(y) and qo y , (see Eqs. (pj|) , (|l7|) ) , can be 
written as follows: 

U+ ^fi W> 2 (J/)) m = (/ <T(r',r)d y ,Q(y')dr'^ </>(y), (29) 

- cf d 2 y Q(y) = p-% J a(r,r')0VHr'. (30) 
After the exclusion of Q(y) from Eq. (pS|), it can be rewritten in the following form: 

(~2^^ + AC/(r) + ^ 2(y) ~/ ^eff(r,rO<W)dr') <t>{y) = ^(y) (31) 

where U e s(r, r') = Jer(r", r)a(r", r')dr"/pcf, U(r) = Jer(r', r)dr', and A = const is defined from the boundary 
condition (m. 

It is easy to see from the structure of j3l| ) that the atom-lattice interaction er(r, r') (when exceeding a certain 
magnitude) can induce an attraction between the gas atoms in the boundary region. Indeed, there can exist such a 
scale of |r — r'| that v5(y — r') — U c s(r, r') < 0. In that case, one would expect essential changes in the spectrum of 



boundary excitations. The study of exciton superfluidity 15 1, EJ hints such a possibility. 



This study reveals a mechanism of the exciton-exciton attraction induced by the lattice effects. The exciton gas falls 
into the soliton-like state ipo, <?0i when the exciton-phonon coupling constant exceeds a certain value, or, equivalently, 
the velocity v exceeds some critical value. The excitation spectrum in an exciton branch has a gap, and the Landau 
critical velocity, calculated for this type of excitations, is given by 

v c ~ K/mL, (32) 

where L denotes the characteristic size of the soliton. We can adapt the similar considerations in our model, replacing 



the exciton-phonon by the gas- wall interactions. The critical velocity can be given by a formula similar to (32) with 
L ~ I, provided I is the only macroscopical length in the theory. 

In this article we did not consider the influence of the long-range van der Waals forces between the walls and gas 
atoms on the stability of superfluid flow. Although the attractive part of these forces originates from interaction 
between the electron shells of the particles] 4lJ , it can be included in our model in the form of a static van der Waals 
potential appearing in the r.h.s. of Eq. (111). Such an external potential being localized near the boundaries can be 
stronger than the effective potential AU(y) in (^lj) and therefore can change the properties of the condensate wave 
function, the boundary excitations, etc.. 

V. CONCLUSIONS 

In conclusion, our simple model manifests one of the possible microscopic mechanisms for dissipation processes in 
the quantum Hamiltonian system of coupled Bose-gas and channel walls. We show that the dissipation can also be 
caused by creation of boundary excitations in this system. Although, in the semiclassical approximation, this process 
is not prohibited at any velocity of the moving repulsive Bose-gas, the higher quantum corrections to the self-energy 
part of the boundary excitations may be essential to obtain the nonzero value of v c in the theory of the dilute Bose- 
gas with gas-wall interaction. On the other hand, more rigorous consideration in the framework of the semiclassical 
approximation should involve solutions of (p|)-(p^) with U ^ 0, where the attractive part of gas- wall interaction is 
taken into account. Such solutions can be still represented in the form (21), (p2[), though equations (|l8|), (|l9|), ( p0| ) 



would become integral. A direct numerical study of the flow dynamics would be also useful. 
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